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1.  Introduction  and  Statement  of  Main  Results. 

Let  M and  N be  random  integers  chosen  uniformly  and  independently 
from  {1,2, . . . ,x) . Throughout  (14, N)  will  denote  the  greatest  common 
divisor  and  [M,N]  the  least  common  multiple.  Cesaro  (1885)  studied  the 
moments  of  (M, W)  and  [M,W].  Theorems  1 and  2 extend  his  work  by  providing 
esrplicit  error  terms.  The  distribution  of  (M,N)  and  [M,W]  is  given  by: 

Theorem  1. 

(1)  P t[M,H]<  tx2  and  (M,U)  =k]  = -|  ~ (kt  (l-log  kt)}  + 0R 

it  k 3 

6 1 log 

(2)  P { (M,N)  = k}  = ~ + 0 (-££-£-) 

it  k 

o k [l/t] 

(3)  P ([M,W]  < tx  } = 1 + -§  * I Id t (1-log  jt)-l)  + 0 (i2^)  . 

« 0=1 

Where  [x]  denotes  the  greatest  integer  less  than  or  equal  to  x. 

Christopher  (1956)  gave  a weaker  form  of  (2). 

(2)  easily  yields  an  estimate  for  the  expected  value  of  (M,N): 

E((M,N))=-|  1 (*>0)  = £ k F _{(M,W)  = k)  = log  x + 0(1). 

x i,d  <x  k < x it 

(2)  does  not  lead  to  an  estimate  for  higher  moments  of  (M,N).  Similarly 
the  form  of  (3)  makes  direct  computation  of  moments  of  [M,N]  unwieldy. 

Using  elementary  arguments  we  will  show: 
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f 


t 

Theorem  2. 

(4)  E { (M, N ) ) = -§  log  x + C + 0 (-&-£) 

x jT 

where  C is  an  explicitly  calculated  constant. 

(5)  for  k > 2,  Ex{(M,N)k}  = - l)  + 0(xk“2  log  x) 

where  £ (z)  is  Riemann' s zeta  function, 

(6)  for  k > 1,  E {[M,N]k}  = ■ x2k  + 0 (x^k  1 log  x)  . 

5(2) (k+1 f 

Section  two  of  this  paper  contains  proofs  while  section  three  contains 
remarks,  further  references  and  an  application  to  the  statistical  problem  of 
reconstructing  the  sample  size  given  a table  of  rounded  percentages. 

2 . Proofs  of  Main  Theorems. 

Throughout  we  use  the  elementary  estimate 

(2.1)  «(x)  - £ cp(k)  - x2  + R(x) 

1 < k < x « 

where  R(x)  = 0 (x  log  x). 

See  for  example  Hardy  and  Wright  (i960)  theorem  360.  Since 

jf  {m,n  < x:  (m,n)  = l)  = 20(x)  + 0(l)  and  (m,n)  = k if  and  only 

if  k|m,  k|n  and  (p  —)  - 1,  we  see  that  Jf=  {m,n  < x:  (m,n)  = k}  = 2$(^-)+0(l). 

This  proves  (2).  To  prove  (l)  and  (3)  we  need  a preparatory  lemma. 
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2. 

Lemma  1.  If  F (t)  = #{m,n  < x:  mn  < t x and  (m,n)  = 1}  , then 

~ X 

yt)  = t(l-log  t)x2  + Ot(x  log  x)  . 

TC 

Proof.  Consider  the  number  of  lattice  points  in  the  region 

2 

Rx(t)  = {m,n  < x:  mn  < tx  }.  It  is  easy  to  see  that  there  are 
2 

t(l-log  t)x  + 0,  (x)  = N (t)  such  points.  Also,  the  pair  (m,n)  e R (t) 

0 X X 

and  (m,n)  = k if  and  only  if  (p  -}  e Rx/fc(t)  (p  f)  = !• 

N (t)  = £ F /,  (t).  The  standard  inversion  formula  says 

x 1 < d < x X/a 

Fx(t)  = £ d(d)N  /d(t)  - ~ t (1-log  t)x2  + 0t(x  log  x)  . 

1 < d < x ' it 

Lemma  1 immediately  implies  that  the  product  of  2 random  integers  is 
independent  of  their  greatest  common  divisor: 

Corollary  1. 

PXUMN]  <tx2j  (M,R)  = k)  = t (1-log  t)  + 0t  k£2jp)  . 


To  prove  (l)  note  that 

P £ [M,N]  < tx2  and  (M,N)  = k}  = P {[[M,N]  < tx2|  (M,N)  = k)  • P {(M,N)w=k) 

X X X 

= P (m  < y x2|'(M,N)  = k)  . P f (M,N)  = k)  . 


Use  of  (2)  and  Corollary  1 completes  the  proof  of  (l).  To  prove  (3)  note 
that 
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[1/t] 


P {[M,N]  < tx2}  = P {( M,U)  > [-]}  + t P {[M,N]  < tx2|  (M,N)=fc)  • P {(M,N)=k} 

X"”  X o n ..  x — X 


L f ■ ■ *4  y ' 

k=l 


Using  (2)  and  Corollary  1 as  before  completes  the  proof  of  theorem  1. 
To  prove  Theorem  2,  write,  for  k > 1, 


(2.2)  £ (m,n)k  = 2 £ £ (m,n)*  £ i 

m,n  <x  l<m<xl<n<m  l<i<x 


.k 


k+1 

= 2 I fk(”)  -lZZ+°(xl 

1 < m < x 


,k 


where  fk(m)  = Y -ru 
i|r 


d cp (— ) . Dirichlets  Hyperbole  argument  (see  eg.  Saffari 

d|m 

(1970))  yields  for  any  t. 


(2-3)  X fk<“>  =.  X > + , X ,»(i)ikf)- y<o®<f> 

l<m<x  1 < 1 < t 1 < 1 < x/t 


where 


i , k+1  , 

ik(t)  = X i -TST  + o(t  > • 


1 < i < t 


Mien  k = 1,  we  proceed  as  follows:  Choose  t - y/x  . The  first  sum 
on  the  right  side  of  (2.3)  is. 


(2.10 


, X +0(|log|)) 

1 < k < /x  n 


= -4-  x2{log  v/x  + 7 + 0(-i)}  + Ofx^2  log  x) 
it2 
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The  second  s\im  in  (2*5 ) is 


(2.5)  £ _q>(H)(i  (f)2  + 0<§)}  =| 

1 < k < /x 


l + 6^) 

1 < k < fX  k 


Now 


where 


£M 


2k+l 


1 < k < -fx  k 1 < k < /x  (k(k+l))‘ 


;r®(k)  + Tu 


2 I 


1 < k < y'x  k(k+l)2  *2 


1 (4  k2  + R(k)}  + l 


$(k) 


1 < k < /x  k2(k+l)2 


+■  -5-  +o  (iSIL-E) 


yx 


2 Z 


Rfr) 


i(k) 


6 r7  k 

2 ^ , , p ' & L>  ' s2  ' & 2 / .2 

it  l<k</x  (k+1)  k=l  k(k-KL)  k=l  k (k+l) 


+ 4 + o(±SS_S) 


Ti 


A 


~r;  log  x + d + 0 


31 


(log  X) 

A 


(2.6)  d = f (<D(k)+2kR(k)  - -§  k(2k+l)),/(k(k+l))2  + -§  (/  + §) 

k=l  Jt  it 


and  7 is  Euler's  constant. 
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Using  this  in  equation  (2.5)  yields  that  the  second  sum  in  (2.3)  is 


(2.7) 


log  X + | X2  + 0(x^2  log  x)  . 

2 it 


The  third  term,  in  (2.3)  is 


(2.8)  i ~ x2  + 0(x5//2  log  x)  . 

it 


Combining  (2.8),  (2.7)  and  (2.4)  in  (2.3)  and  using  this  in  (2.2)  yields: 


X (m,n)  = x2  log  x + (d  + (y  + |)  - jy)x2  + 0(x3;/2  log  x)  , 
m,n  < x « it 

■where  d is  defined  in  (2.6). 

Then  k > 2,  the  best  choice  of  t in  (2.3)  is  t = 1.  A calculation 
very  similar  to  the  case  of  k = 1 leads  to  (3). 

We  now  prove  (6).  Consider  the  sum 


(2.9) 


I [i,Uk-2  l l [i,3)k  + 0(xk+1) 


i,j  < x 


i < x j < i 


= 2 


l % t (¥>  +o(xk+1) 

i < x d|i  3 < i 


= 2 l lk  I f .(£)  ♦ 0(xk+1)  - 2 £dk  l i\(3)*0  (xk+1) 
i < x d|i  d=l  3 < x/d 


Where 
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yn)  - 


E 3 

3 < n 

(J  j»n)=l 


,k 


We  may  derive  another  expression  for  (n)  by  considering  the  sum 


(2.10) 


n , k+1 

X,  1 'Vdt Vn>  - 

x=l 


n 


k r>  k 


d | n dk 


Dividing  (2.10)  by  n and  inverting  yields 


fk^  1 y 

k ~ k+1  V ^d^  ? ^ ^ k 

n d n d n d 


or 


Vn>  = ra:  *<n>  * J Vd)  * ^ + E(”)  • 

d|n 

When  we  substitute  this  expression  for  f^(j)  in  (2.9)  we  must  evaluate: 

Sl(y)  = . § ^(j)  = . £ 3 ?.  ^d^db 

3 <y  j < y a|j 

= E ^(i)i2k  £ \(d)dk  . 

i < y d < y/i 

Wow  R^(d)  Is  a polynomial  in  d of  degree  k.  Thus, 


I spy)  I < 


I i2k(f)2k+1 

• - -L 

x < y 


n t 2k+l  . 

0 (y  log  y) 
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We  must  also  evaluate 


S2(y)=k^I  E A(D)=£^£2k  l +o(  l D2k"23>(j))+$(y)y2*} 

j < y d < y d <y 

OVJ-p 

6 k y , 3 l . 2k+2  , „ , 2k+l  , 

-7  Tt*rj  ~W*T)  7 (51) y 0(y  logy) 

6 A k u 2k+2  , „ , 2k+l  , \ 3 1 2k+2  . „ , 2k+l  , 

= ~2 — : (p  ~ + °(y  lQg  y)  = -5 o y +°(y  log  y) 

* k+l  2 d*  Z n2  (k+1)2 


Substituting  in  the  right  side  of  (2.9)  we  have 


l [1,3]*-*  I a*  (Si(|)-«-s2(|)}  +o(xk+1) 

i , d < x d=l 


6 1 2k+2  £ 1 , w 2k+l  . 

= “p  ~o  x 2*  -T-fo+0^x  log  x) 

It  (k+1)  d=l  d 


= g(k+2)  x 


2k+2 


f(2)  (k+1/ 


, n / 2k+l  _ \ 

+ 0 (x  log  x ) . 


Miscellaneous  Remarks. 

1.  If  . . . ,M^.  are  random  integers  chosen  uniformly  at  random 

then  results  stated  in  Christopher  (1956)  (see  also  Cohen  (i960),  Herzog 
and  Stewart  (1971  )>  and  Neymann  (1972))  imply  that 


(3.1) 


F*f<vv-"’V^>-db  i + 0(r^i) 


x 3 


We  have  not  tried  to  extend  theorems  1 and  2 to  the  k-dimensional  case. 


8 


3.1  has  an  application  to  a problem  in  applied  statistics.  Suppose 
a population  of  n individuals  is  distributed  into  k categories  with 


n individuals  in  category  i.  Often  only  the  proportions  p^  = n^/n 
are  reported.  A method  for  estimating  n given  p^,  1 < i < k is 
described  in  Wallis  and  Roberts  (1956),  pgs.  18U-I89.  Briefly,  let 


m = min 


I 

i=l 


where  the  minimum  is  taken  over  all  k tuples 


(b1,b0, . . . ,bk),  with  b^  € {0,  +1,  +2,...  } not  all  V equal  zero. 

An  estimate  for  n is  [l/m].  This  method  works  if  the  p\  are  reported 
with  enough  precision  and  the  n.  are  relatively  prime  for  then  the 
Euclidean  algorithm  implies  there  are  integers  (b^K_^  such  that 

XI 

b.n.  = 1.  These  b.  give  the  minimum  m = — . If  it  is  reasonable 
11  1 n 

to  approximate  the  n^  as  random  integers  then  (3-l)  implies  that 
Prob  ( (n^,n2, . . . ,nk)=l)  = y^-y  as  expected,  as  k increases  this 

probability  goes  to  1.  For  example,  = .96^,  ^yyry  = *992> 


Wi 


= .998*  This  suggests  the  method  has  a good  chance  of  working  with  a 


small  number  of  categories.  Wallace  and  Roberts  (1956)  give  several  examples 
and  further  details  about  practical  implementation. 


2.  The  best  result  we  know  for  R(x)  defined  in  (2.1)  is  due  to 
Saltykov  (i960).  He  shows  that 

R(x)  = 0 (x (log  x)2/^(loglog  x)1+€)  . 

Use  of  this  throughout  leads  to  a slight  improvement  in  the  bounds  of 
theorems  1 and  2. 

3*  The  functions  (M,W)  and  [M,N]  are  both  multiplicative  in  the 
sense  of  Delange  (1969,1970).  It  would  be  of  interest  to  derive  results 
similar  to  Theorems  1 and  2 for  more  general  multiplicative  functions. 

Q 
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